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Abstract: The major contribution to the anisotropy of the temperature of the Cosmic 
Microwave Background (CMB) radiation is believed to come from the interaction of lin- 
ear density perturbations with the radiation previous to the decoupling time. Assuming 
a standard thermal history for the gas after recombination, only the gravitational field 
produced by the linear density perturbations present on a ft ^ 1 universe can generate 
anisotropics at low z (these anisotropies would manifest on large angular scales). How- 
ever, secondary anisotropies are inevitably produced during the nonlinear evolution of 
matter at late times even in a universe with a standard thermal history. Two effects 
associated to this nonlinear phase can give rise to new anisotropies: the time-varying 
gravitational potential of nonlinear structures (Rees-Sciama RS effect) and the inverse 
Compton scattering of the microwave photons with hot electrons in clusters of galaxies 
(Sunyaev-Zeldovich SZ effect). These two effects can produce distinct imprints on the 
CMB temperature anisotropy. We discuss the amplitude of the anisotropies expected 
and the relevant angular scales in different cosmological scenarios. Future sensitive ex- 
periments will be able to probe the CMB anisotropies beyong the first order primary 
contribution. 



1 Introduction 



Cosmic Microwave Background (CMB) temperature anisotropies carry informa- 
tion on the large scale distribution of matter via the gravitational imprint of 
linear density fluctuations on large angular scales ( ^ [2fi 1 / 2 )°), the so called 
Sachs- Wolfe effect for a SI = 1 universe (Sachs and Wolfe 1967) or the general- 
ized Sachs- Wolfe effect for a Q ^ 1 universe or in the presence of a cosmological 
constant A (Anile and Motta 1967, Wilson 1983, Gouda et al. 1991). In the 
first case the gravitational potential is static and the temperature fluctuations 
are produced by the potential fluctuations at recombination whereas in the other 
cases they are also produced by an integrated effect of the time- varying potential 
along the photon path from recombination to the observer. 
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On smaller angular scales the major contribution to the CMB anisotropy 
comes from photon fluctuations and the Doppler effect due to the velocity of the 
last scatterers at recombination. These primary anisotropies can be erased if the 
gas in the universe is reionized at a sufficiently early time (z ^ 100). However, 
in this case new secondary anisotropies are generated by the gas motion which 
produces a Doppler effect and by the coupling of velocity and density fluctuations 
(Vishniac f 987). All the mentioned effects characterize the shape of the radiation 
power spectrum which contains information on the cosmological parameters Q, 
A, the baryonic content the Hubble parameter, the initial matter density 
distribution and the thermal history of the universe (see Hu in this volume for 
a general review on the CMB anisotropies). 

CMB temperature anisotropies, detected for the first time by Smoot et al. 
(1992) with the COBE-DMR experiment, are believed to be dominated by the 
interaction of matter density perturbations and radiation to first order in per- 
turbation theory. In the standard scenario these primary anisotropies are pro- 
duced before neutral hydrogen formed at recombination (z w 1000) and they are 
preserved until today. However, secondary anisotropies are inevitably produced 
during the nonlinear evolution of matter at late times. Even in a flat universe 
with a standard thermal history (i.e. with no reionization of the gas content af- 
ter recombination), the time- varying gravitational potential due to the evolving 
nonlinear structures generates new temperature fluctuations (Rees-Sciama RS 
effect, Rees and Sciama 1968). Their amplitude is estimated to be ~ 10 -6 and 
the signatures imprinted in the radiation power spectrum may distinguish them 
from the primary anisotropies. We will discuss the possible detectability of these 
secondary fluctuations with future very sensitive experiments. 

The Compton scattering of microwave photons by hot electrons in the intra- 
cluster gas also produces new anisotropies at very late times after the formation 
of rich clusters of galaxies (Sunyaev and Zcldovich 1972). This Sunyacv-Zcldovich 
effect (SZ) produces spectral distortions which translate into anisotropies when 
the angular distribution of temperatures in the sky is measured by an experi- 
ment at a given frequency band. The peculiar velocity of clusters also generates 
anisotropy via Doppler effect but the temperature amplitude of this kinematic 
Sunyaev-Zeldovich effect (KSZ) is expected to be at least an order of magnitude 
below the thermal SZ one in the Raileygh- Jeans region of the spectrum. The SZ 
effect is several orders of magnitude greater than the RS effect for a rich cluster. 
However, since the SZ effect is restricted to the central regions of rich clusters 
with a high concentration of hot gas whereas the RS one is a much more general 
effect the statistical average of many structures in the sky may give comparable 
amplitudes of the two effects or may even be dominant for the later effect. 



^f- Beyond Linear Theory 



2 The Rees-Sciama Effect 



Rees and Sciama (1968) first noticed that the time-varying gravitational poten- 
tial of evolving nonlinear structures will produce secondary anisotropies in the 
temperature of the CMB. That seminal paper, where an incorrect Newtonian 
approach was used to study the anisotropy generated by a compensated cluster, 
motivated many works which used a General Rclativistic treatment to model the 
nonlinear evolution of spherical over-densities and under-densities in a uniform 
background. A simple model used was that of the "Swiss-cheese" where it is 
assumed a central spherical over-density surrounded by a compensated under- 
density (Notale 1984, Martinez-Gonzalez and Sanz 1990). Also a similar model 
consisting of a central spherical under-dense region compensated by a thin shell 
was used for cosmic voids (Thompson and Vishniac 1987, Dyer and Ip 1988, 
Martinez-Gonzalez and Sanz 1990). A more sophisticated spherical model based 
on the Tolman solution of the Einstein equations allows a general density profile 
(Panek 1992, Arnau et al. 1993, Fang and Wu 1993, Fullana et al. 1996). 

A useful expression to calculate the RS effect produced by non-spherical 
structures and general realistic matter distributions can be derived within the 
potential approximation (Martinez-Gonzalez et al. 1990, see below). Chodor- 
owski (1992, 1994) and Atrio-Barandcla and Kashlinski (1992) used that ex- 
pression to calculate the effect of pancake-like collapse of ellipsoids. The results 
of all these works for the most prominent structures seen in our local universe 
can be summarized as follows. For the large Bootes void of size « 60/i _1 Mpc 
located at a distance of ss 150 hr 1 Mpc and with a density of galaxies w 25% 
of the background density (Kirshner et al. 1981, de Lapparent et al. 1986), it 
is estimated AT jT w a few x 10~ 7 with small variation for different models 
(Martmez-Gonzalez and Sanz 1990, Panek 1992, Arnau et al. 1993). For the 
Great Attractor located at a distance of « 60 h~ x Mpc and a density contrast 
of w 1 inside a radius of « 40/i _1 Mpc (Lynden-Bell et al. 1987) the estimated 
RS effect is « a few x 10~ 6 and again the results are fairly insensitive to the 
model used (Martmez-Gonzalez et al. 1990, Panek 1992, Arnau et al. 1994). In 
the case of the large concentration of galaxies known as the Great Wall which 
forms a flattened structure of dimensions « 5 x 60 x 170 hr x Mpc at a distance 
of some 75/i _1 Mpc (Geller and Huchra 1989) the effect is also a few x 10~ 6 
(Atrio-Barandela and Kashlinski 1992, Chodorowski 1994). 

If similar structures are assumed to be at different distances from the observer 
in low density universes then it is possible to increase the RS effect in almost 
an order of magnitude. For distances z ~ 1 — 10 and universes with Q « 0.2 
great attractor-like structures can produce an effect of the order of 10 -5 (Arnau 
et al. 1994) and Bootes-like voids of a few x 10~ 6 (Fullana et al. 1996), with 
angular scales of a few degrees. Thus, only in open universes the presence of such 
large inhomogeneities might leave an observable imprint in the CMB anisotropy 
maps. 
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Below we will describe a useful expression for the RS effect in realistic situ- 
ations and use it to calculate the statistical effect produced by a distribution of 
matter. 

2.1 The Potential Approximation 

Martinez-Gonzalez et al. (1990), based on the potential approximation to gen- 
eral relativity, derived a useful expression for the temperature fluctuation gen- 
erated when microwave photons cross a nonlinear density perturbation A(t, x). 
If ip(t, x) is the corresponding gravitational potential then the RS effect can be 
simply written as 

( AT / T ) secondary^ 2 £ dt ^ X ) > V*<P = Ma^Afrx) , (1) 

where a(t) is the scale factor normalized to the present time (a — 1) and our 
units are c = 8nG = 1 and the Hubble length at present is d a = 2/H~ 1 . 
Thus, the Rees-Sciama effect is basically the work performed by a photon trav- 
elling through the time-dependent gravitational potential ip(t, x) from recom- 
bination to the observer. The previous equation was initially derived for a flat 
universe, however, since the scales relevant for nonlinear structures are always 
much smaller than the curvature scale it also applies to open universes. The 
integral must be performed along the geodesic associated to the corresponding 
background 

x(t, n) = A(a)n , A(a) = 1 - (1 _^ /a + , L n) -V> ^ 

An order of magnitude estimate of the effect can easily be obtained from 
equation (1). For a photon crossing a cluster the anisotropy is 




where <p is the potential of the cluster, t c and td are the crossing time and the 
dynamical time and v the velocity of collapse. For a cluster with v = 1000 km/s 
the effect is AT jT ~ 10~ 7 . A statistical distribution of lumps can therefore 
increase the effect to values of cosmological interest. 

2.2 Second Order Perturbation Theory 

We will now estimate the RS effect in an open or flat background for a cold dark 
matter (CDM) cosmogony. This effect has been recently calculated by Sanz et 
al. (1996) within the second order perturbation theory. We will follow that paper 
below. 

For vanishing pressure, by perturbing the Einstein field equations and consid- 
ering perturbations in density up to the second order we can obtain the following 
expression for the density fluctuations 
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A(t,x) = D6 + D 2 [p 2 + W6-V^+^ ij ^] , v 2 e = 5 (4) 

where S(x) = 5 r (l + z r ). Equation (4) is exact within the 2nd order perturbation 
theory for a fiat background (see Peebles 1980) and it is a good approximation 
for open universes with J? Z, 0.1 as has been shown by Bouchet et al. (1993) and 
Catelan et al. (1995). By using equations (2,3) in the expression for the second 
order effect (eq. 1), we obtain the following result 

AT f° 

( T 1 secondary d\W{X)(f)(x \fl) , 

J Ms 

W(\)=24f2D 2 (2f- 1) 1 ^^ (5) 

where the new potential related to the second order desity perturbation <$ 2 sat- 
isfies the following Poisson equation 

V 2 #r)^ 2 ^ 2 + V^V£+^e 4 -' • (6) 

In the equation above D is the growing mode of the perturbations normalized to 
1 at the present time. For a flat universe: D — a, whereas for an open universe 
(Peebles, 1980): 

g{x ) U din a 

g(x) = l + (-)\l + (l + -) 1 / 2 \n((l+x)^ 2 -x^ 2 )] . (7) 
X L x \ /J 

It is common to expand the temperature fluctuations in terms of spheri- 
cal harmonics Yi m , AT(n) — Sm=-i a lmYi m (p>)- The temperature corre- 
lation function is given in terms of the radiation power spectrum as C(<9) = 
l/(47r) ^2(2l + l)CiPi(cos(a)), with Ci =< af m >. In the second order perturba- 
tion theory the multipole component Ci is (we shall not consider the monopole 
and dipolc in the calculations of the temperature anisotropy below) : 

Q = 1J dkk- 2 P {2) (k)R?(k) ,R l {k) = £" dXW(X)ji{kp) . (8) 

Here, p = A/[l — (1 — J?) A 2 ], ji is the Bessel function of fractional order and 
Xis is the distance from the observer to the last scattering surface. This equa- 
tion generalizes the result obtained by Martinez-Gonzalez et al. (1992) to open 
universes. The function P( 2 ) (k) is the power spectrum associated with the 2nd 
order density perturbation <5 2 and is related to the power spectrum P^ of the time 
derivative of the potential £ by P^ = prP( 2 )(fc). The second order perturbation 
power spectrum P^ 2 ) (k) is given in terms of the first order power spectrum P(k) 
by the equation (Goroff et al. 1986, Suto and Sasaki, 1991) 
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k 3 f 00 f 1 
P(2)(fc) = 98(2^ J drP{kr) J_, dxP ^ r2 + 1 - 2rx ) V2 ) 

/3r + 7x - 10ra; 2 \ 2 

V r 2 + l-2rx ) ' ^ 
In the limit of small k, there is a cancellation of the three terms contributing to 
62, implying that P( 2 ) (k) has very little power on large scales and in this regime 
goes like P( 2 )(fc) oc k A independently of the primordial power spectrum. 

Equation (8) for the multipoles Ci can be further simplified by realizing 
that the integral to calculate Ri{k) is the product of a slowly changing function 
of A, the window W(X), and a spherical Bessel function (Seljak 1996). Thus, 
fixing the value of the argument of W and applying the large I approximation 
Jo Ji( x ') — V^ 77 '/2Z)5(x — I), S being the step function, we obtain 

c t = ^ £' s dxw 2 (x)p 2 [i(i-(i-n)x 2 )/\] . (10) 

This approximation works very well and in the worst situation for the 
quadrupole moment the error is ^ 30% for all Q values. 

In relation to a possible secondary contribution coming from 3rd order density 
perturbations for open models (in the case of flat models this does not exist 
because the 1st order gravitational potential is static), we have the following 
comment: the coupling of the 1st order gravitational potential tp^ oc — with 
the 3rd order potential ip^ oc — gives a kernel for the integrated gravitational 
effect proportional to -§r(/ — 1)(3/ — 1), whereas the coupling of the 2nd order 
gravitational potential ip^ oc ^- with itself gives a kernel for the integrated 

gravitational effect proportional to ^r(2/ — l) 2 . This second function is always 
greater than the first one. Moreover, for quasi-flat models we expect a negligible 
contribution from the coupling of lst-3rd order perturbations because / « 1, 
whereas for low-.!? models the integrated gravitational effect due to 2nd-2nd 
order perturbations is produced at high-z (75% of the final effect is produced in 
the interval [10, 10 3 ] for fl = 0.1, see next section) where / w 1, and so there 
is practically no lst-3rd order contribution. At smaller z some contribution due 
to the lst-3rd coupling is produced but it is estimated to be always bounded by 
that due to the 2nd-2nd coupling at low-z, and this is a small fraction of the 
final contribution. 



2.3 Results 

The equations derived in the previous section have been applied to calculate 
the predicted amplitudes of the multipole components C; in open universes. 
We assume an open or flat CDM model with a Harrison-Zeldovich primordial 
spectrum of matter density fluctuations P(k) = Ak and J?b = 0.05, H = 50km 
s _1 Mpc -1 . In figure 1 we display the multipoles C; for fi = 1,0.3,0.1 due 
to primary anisotropics (upper curves) and due to 2nd order anisotropics (lower 
curves). Primary anisotropics are normalized to the 2-year COBE-DMR maps as 



4f- Beyond Linear Theory 7 



given by the analysis of Cayon et al. (1996) for the HZ spectrum. Since secondary 
anisotropics are generated by the relatively small scale structure, ^ lOOMpc, 
it is more appropriate to use the a w — 1 normalization, i.e. the rms density 
fluctuation is unity at 16 Mpc. The secondary radiation power spectrum peaks 
at I « 250 for all fl values, contrary to the linear one where the multipole order, 
Z, of the peak increases for low fl values. For the secondary contribution, the shift 
of the maximum of the 2nd order power spectrum P( 2 ) (k) towards large scales 
(small k) for low fl values is balanced by the opposite shift of the maximum of 
the function Ri(k) in equation (8). 



Fig. 1. Radiation power spectrum for the primary (upper curves) and secondary (lower 
curves) contributions. Solid, dashed and dotted lines correspond to O — 1,0.3,0.1, 
respectively. The secondary contribution is normalized to (Tie = 1, except for the solid 
upper curve which represents the flat case with the COBE-DMR normalization. 



The generation of the secondary radiation power spectrum with redshift is 
shown in figures (2a, b) for fl = 1,0.1 respectively. For fl = 1 more than 90% 
of the anisotropy is produced at z ^ 10 (Martinez-Gonzalez et al. 1992). In the 
case of low fl models the anisotropy is produced at a relatively higher redshift, 
for fl = 0.1 about 80% of the effect is generated at a redshift z < 30. In any 
case, reionization of the matter in the universe cannot substantially erase the 
secondary anisotropy at such low redshifts. 

Nonlinear scales do not enter in the second-order calculation (eq. 8) because 
for any fl value the maximum of Ri(k) is at k w Ifl and the scales that are 
contributing to the multipole I have k I. Thus, for the multipolcs C; up to 
I = 1000 only scales > 12Mpc contribute to the second-order effect. This result 



8 Enrique Martinez- Gonzalez 



Fig. 2a. Generation of the secondary radiation power spectrum with redshift for a flat 
universe. Solid, dashed and dotted lines correspond to the anisotropy generated from 
redshifts z = 100, 10, 1 to the present, respectively. 



Fig. 2b. The same as figure (2a) but for an open universe with = 0.1. 
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agrees with the recent calculations by Seljak (1996) considering nonlinear scales, 
up to I ~ 1000 the RS effect is dominated by the second-order calculations. 
Tului and Laguna (1995) and Tului et al. (1996) have computed the RS effect 
by propagating a bundle of photons through an evolving N-body simulation 
from the decoupling to the present time. However, one must be careful when 
calculating the nonlinear effect because at scales £ lOMpc the use of eq. (1) to 
compute the RS effect may not be correct. The reason is that at those small 
scales higher orders including spacial derivatives of the gravitational potential, 
which have been neglected in the perturbative expansion of the Einstein field 
equations to derive eq. (1), may contribute significantly. 

The results mentioned above are for the two-point temperature correlation 
function where clearly the primary contribution dominates over the secondary 
RS one. However, the linear primary effect gives a null contribution to the three- 
point correlation since the primordial matter perturbations generated during the 
inflationary phase are Gaussian. Recently, Munshi et al. (1995) and Mollerach et 
al. (1995) have computed the leading contribution to the three-point correlation 
function which involves the second-order RS effect. Unfortunately the predicted 
signal is below the unavoidable cosmic variance uncertainty. 

3 The Sunyaev-Zeldovich Effect 

Clusters of galaxies not only induce secondary anisotropics in the CMB by means 
of gravitational effects. As first noticed by Sunyaev and Zeldovich (1970,1972) 
the scattering of microwave photons by hot electrons in the intracluster gas 
produces spectral distortions in the blackbody spectrum, which is known as the 
thermal SZ effect. Additionally, the peculiar motion of the clusters also produces 
anisotropy in the CMB via Doppler effect, known as the kinematic SZ effect. In 
subsections 3.1 and 3.2 we will study in detail these two effects (see also the 
recent review by Rephaeli 1995). 

The thermal SZ effect has been measured in a number of rich Abell clusters 
using three distinct techniques: single-dish radiometry, bolometric observations 
and interferometry (see Birkinshaw 1993). Some of the first convincing detections 
were those of the clusters A665, A2218 and 0016+16 using the single-dish OVRO 
radio telescope at 20GHz (Birkinshaw et al. 1984). More recently Herbig et al. 
(1995) have measured the effect in the nearby Coma cluster at a high significant 
level with the OVRO 5.5m telescope at 32GHz, deriving a maximum Rayleigh- 
Jeans decrement of AT RJ = —505 ± 92/xK. The first detection of the SZ effect 
in the millimetric region has been obtained with the bolometric observations of 
the SUZIE experiment on A2163 at 2.2mm, close to the wavelength at which 
the SZ decrement is maximum (Wilbanks et al. 1994). With these bolometric 
observations it is expected to measure the positive effect at wavelengths < 1.4mm 
(as we will see below the thermal SZ effect changes sign at that wavelength). The 
first image of the SZ effect was obtained for A2218 with the Ryle interferometer 
at 15Ghz (Jones et al. 1993). More recently, images of the rich clusters A773 
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and 0016+16 have been taken with the Ryle and OVRO arrays (Grainge et al. 
1993, Saunders 1995, Carlstrom et al. 1996) and several other clusters are being 
imaged at present (very recently, an image of the SZ decrement towards A1413 
has been taken by Grainge et al. 1996 with the Ryle Telescope). Interferometry 
makes possible to directly compare the SZ image with the X-ray image of a 
cluster and thus study the structure of the intracluster gas. 



3.1 Thermal SZ Effect 

The inverse Compton scattering of the isotropic CMB radiation (hv « m e c 2 ) 
by a hot (T e >> T), nonrelativistic (kT e « m e c 2 ) Maxwellian electron gas 
is quantitatively described by a simplified version of the Kompaneets (1957) 
equation for the rate of change of the photon occupation number n (Sunyaev 
and Zcldovich 1980) 

dn 1 d 4 dn 
dy x 2 dx dx ' 

kT e hv 
^dr , x = — 

m e c 2 kT 

where the nondimensional frequency x has been defined, the optical depth due 
to Thompson scattering is dr = n e <TTcdt with ot being the Thomson cross 
section and y is the Comptonization parameter. Since for the richest clusters 
the radiation is only weakly scattered, y « 1, and so the deviations from the 
Planckian spectrum are small, then it is easy to solve equation (11) by inserting 
in the right-hand side of this equation the Planck function for the occupation 
number np = l/(e x — 1) (Zcldovich and Sunyaev 1969). The solution for the 
change of spectral intensity I = 2hv 3 n/c 2 is 



2(kT) 3 x A e x 

\2 ■ 



Al 



X 



x coth 

2 



(he) 2 (e x -lf 

The corresponding change in the thermodinamic temperature is 



(12) 



AT = yT[x coth^ -4] . (13) 

For the Raileygh- Jeans region of the spectrum, x « 1, equations (12,13) 
take the following simple form: 

AIrj ~ "Tkp" x2y > at rj-- 2 v t - ( 14 ) 

Previous solutions are obtained from the Kompaneets equation which is based 
on a diffusion approach and thus is not in principle suitable for Compton scat- 
tering in clusters where the optical depth r << 1 and consequently most of 
the microwave photons are not scattered even once. Another limitation is the 
nonrelativistic treatment used which affects to the calculations for the richest 
clusters with temperatures kT e <~ lOkeV and produces significant errors in the 
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Wien region (Fabbri 1981, Raphacli 1995). Even so, except for the extreme cases 
in the distribution of clusters, solutions (12,13) provide an adequate description 
of the SZ effect in the whole spectral range of interest (for the richest clusters 
the solutions are still accurate in the Raileygh- Jeans region). 

An interesting property of equations (12,13) is that the spectral change is 
zero at A « 1.4mm (y w 220GHz) being negative and positive above and below 
that wavelength, respectively. This characteristic change in the sign of the SZ 
effect is still waiting to be observed (a tentative detection of the sign change has 
been recently claimed by Andreani et al. 1996 for the clusters A2744 and S1077 
observed at 1.2 and 2 mm with the 15m SEST antenna in Chile). 

We can now estimate the amplitude of the effect for a cluster of galaxies. 
Considering the nearby, well-studied Coma cluster, we know from X-ray data 
that the properties of the cluster atmosphere consist of a central density of 
n e w 2.9 x lCT 3 cm- 3 , a radius r c w 10'.5 (« 0.43Mpc for H = 50Km s" 1 
Mpc^ 1 ) and a temperature T e w 9.1kcV (Hughes et al. 1988, Briel et al. 1992). 
From these values we can easily derive an approximate optical depth of r « 
n e aT^r c « 5.1 x 1CP 3 and a maximum decrement in the Raileygh- Jeans region 
of AT w 500^K, very close to the value derived from the observations with the 
OVRO 5.5m telescope (Herbig et al. 1995). 

Combining the SZ measurement with the parameters of the model atmo- 
sphere of the cluster derived from the X-ray data one can determine the angular 
diameter distance to a cluster and thus deduce the Hubble constant H D with- 
out relying on any cosmic distance ladder. Estimates of H Q have already been 
made for several clusters of galaxies providing an accuracy for H Q similar to other 
methods (see for instance the recent estimate from the SZ measurement in Coma 
by Herbig et al. 1995). The major uncertainty of this method originates from 
systematic effects related to the X-ray model, departures from spherical symme- 
try and possible existence of dumpiness in the gas distribution (Birkinshaw et 
al. 1991). 



3.2 Kinematic SZ Effect 

The peculiar velocity of a cluster also modifies the intensity of the CMB via the 
Doppler effect. Because of its nature, this kinematic SZ effect produces fluctua- 
tions in the CMB temperature which are independent of frequency, in contrast 
to the thermal SZ. The change in temperature is therefore given by 

AT k = -tT^ (15) 

c 

where v r is the radial velocity of the cluster with positive sign for recession. 
The change in the intensity AI k can be easily obtained from the temperature 
fluctuation using the Planck spectrum 

2(kT) 3 x 4 e x v r 
= ~ (he) 2 { e*-l? T -- < 16 > 
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For a rich cluster like Coma and assuming a peculiar velocity of 500km s^ 1 the 
expected change in the temperature is only of ATk ~ 20^iK, much smaller than 
the thermal effect. It is interesting to notice that the peculiar velocity of a cluster 
can be determined by measuring its thermal and kinematic effects and knowing 
the spatial distribution of the gas temperature from X-ray spectral data. Future 
very sensitive multifrequency ground-based, balloon and satellite experiments 
surveying large areas of the sky at a resolution of £ 10 arcmin (see e.g. the 
proposed satellite COBRAS/SAMBA to ESA M3, Presentation of Assesment 
Study Results 1994) will be able to measure the thermal and kinematic SZ 
effects in ~ 1000 clusters. This information combined with X-ray spectral data 
of the same clusters will allow a reliable measurement of the bulk motion in a 
volume of ~ 1 Gpc (Haehnelt and Tegmark 1995). 

3.3 Contribution to the CMB fluctuations 

In the previous subsections we have studied the SZ effect expected when the 
microwave photons cross the intracluster gas of hot electrons in a rich cluster of 
galaxies. Below we will discuss the contribution of the emsemble of clusters to 
the CMB temperature maps. 

Considerable work has been done to estimate the rms temperature fluctua- 
tions originated by the SZ effect. The usual approach is to consider an evolution 
for the cluster mass function that corresponds to a previously assumed theo- 
retical model of density perturbations (Cole and Kaiser 1989; Markevitch et al. 
1991, 1992; Makino and Suto 1993; Bartlctt and Silk 1994; Colafrancesco et al. 
1994; De Luca et al. 1995). The cluster mass function can be represented by the 
Press-Schechter formula (Press and Schechtcr 1974) which gives the comoving 
number density of collapsed objects per mass interval at a given redshift. A col- 
lapsed object is formed when its density contrast reaches the value 5 C — 1.68 
as given by the spherical collapse model. Assuming that the cluster is virializcd 
to a given virial radius R v the temperature of the gas can be related to the 
mass of the cluster. The electron density profile is usually parametrized with 
a /3-model with the exponent (3 « 0.75 from observations. Typical rms values 
of the temperature fluctuations obtained with this formalism are of the order 
AT/T ~ 10~ 6 . 

A different approach has been considered by Ccballos and Barcons (1994) 
who used an empirically based model for the mass function of the intracluster 
mass. They use a parametrization of the X-ray luminosity function and its neg- 
ative evolution (Edge et al. 1990) which supports the idea that the number of 
luminous clusters decreases with z (Giogia et al. 1990; Henry et al. 1992) and 
therefore the amount of hot gas available for the inverse Compton scattering of 
the CMB photons is consequently limitted. The result is a negligible SZ contri- 
bution to the rms temperature fluctuations of AT jT ^ 10~ 7 . The contribution 
is dominated by the hottest clusters which are easily identified and thus its SZ 
signal can be removed from the temperature map. 

The evidence for fewer X-ray luminous clusters at z > 0.2 than locally shown 
by the Einstein Extended Medium Sensitivity Survey (EMSS) should be con- 
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firmed using a complete, X-ray selected sample. The present ongoing cluster 
surveys based on ROSAT data will measure the X-ray luminosity function at 
redshifts > 0.2 and thus will determine the amount of hot gas at higher redshifts. 
Evidence of the marked decline in the volume density of luminous clusters with 
redshift has been recently presented with the X-ray RIXOS survey up to z w 0.6 
by Castander et al. (1996) (this claim is, however, based on a sample of only 13 
clusters). 



4 Summary 



The RS and SZ effects, associated with the late, nonlinear phase of structure 
formation, can give rise to secondary anisotropies. The RS effect is produced by 
the time- varying gravitational potential of the evolving matter density perturba- 
tions whereas the SZ one is due to the inverse Compton scattering of microwave 
photons by hot electrons in the intracluster gas. These two effects produce dis- 
tinct signatures on the CMB temperature anisotropy: contrary to the RS effect 
the SZ one is frequency dependent. 

The secondary anisotropy generated by the SZ effect is probably dominated 
by the more luminous clusters and is generated at low redshifts z ^ 1, if the neg- 
ative evolution in the X-ray luminosity function of clusters with redshift shown 
by the EMSS is confirmed. The amplitude of the signal is very sensitive to the 
amount of hot gas available at moderate and high reshift, z ^ 0.5. Thus, predic- 
tions based on theoretical models like the CDM one, which imply an increase in 
the amplitude of the X-ray luminosity function with redshift, give anisotropies 
(AT/T) SZ ~ 10~ 6 whereas calculations which take into account the observed 
decline of luminous clusters with z suggest much smaller values. 

The amplitude of the RS effect for realistic models of structure formation is 
{AT/T)rs ~ 10~ 6 , being higher for universes with high fi. The signal is very 
sensitive to the normalization of the matter power spectrum being directly pro- 
portional to its amplitude. Early reionization of the matter in the universe, which 
can strongly affect the primary anisotropies, would not appreciably change the 
secondary ones. Future very sensitive satellite experiments might be able to test 
the gravitational instability theory beyond the linear theory. 
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